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Abstract, T decays into hadrons have a twofold interest: On the one hand, they are a clean 
m : environment for studying the hadronization of the left-handed current of QCD, while, on the 

■^j- ' other side, provide relevant dynamical information of the resonances that mediate these processes. 

Within an effective field theory-like framework, namely Resonance Chiral Theory, we analyse 
the decays ot the z into KKn modes and compare the results with CLEO and BaBar data. In this 
way, we provide bounds on the couplings entering our Lagrangian and predict the corresponding 
spectral functions. As a main result -and contrary to the bulk of theoretical studies and experimental 
analyses- we find vector current dominance on these decays. 
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1. INTRODUCTION 



> 
m 

Although it is perfectly known that Quantum Chromodynamics |T] (QCD) is the quan- 
tum field theory that describes the strong interaction, we do not know how to solve 
it when it becomes non-perturbative. That happens, roughly speaking, below 2 GeV. 
Ideally, one would like to complement partonic QCD with a dual theory written in 
terms of the relevant degrees of freedom in the low and intermediate energy domain 
(mesons and baryons) that could be treated perturbatively at large distances, easening 
the computations and ensuring convergence. Unfortunately, it is not yet known how to 
5-j \ accomplish this task. 

T decays into hadrons allow to study the hadronic properties of vector and axial- vector 
QCD currents and, accordingly, to find the dynamics generated by the exchanged reso- 
nances. At very low energies (E <C M p ), Chiral Perturbation Theory (#PT) |UH|] is the 
effective field theory of QCD, but only enables to explain hadronic tau decays in a tiny 
window of the available phase space |4|]. However, we cannot neglect the importance of 
the resonance states exchanged in the process. Applying vector meson dominance 
the most relevant resonances are the p(770), the »(782), the iT(892) and the ai(1260). 
To explore them, one should include these degrees of freedom explicitly into the theory. 
There have been several strategies to tackle the problem: the best-known one is that 
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of employing a parameterization using a combination of Breit-Wigner's [HI, @]. Even 
though it is clear that polology demands 'something of the kind', it is not clear how to 
go beyond and formalize its link with QCD, if it exists. In fact, the analysis of the decay 
T — ► 37TV T undertaken in Refs. [0, M\ showed that this procedure violates the chiral 
symmetry of massless QCD at 0(p 4 )- 

Several experiments have been collecting good quality data on x — > KK%v x decays. 
In particular, CLEO and BaBar have published the rates for these channels. 

It is known [@] that the decay x~ — > K + K K v T is not well described by a model 
generalizing the one given in Ref. |6{]. As a result of this, the CLEO collaboration 
proposed ft 1 OTI a modification of it, able to fit the data at the price of violating the 
Wess-Zumino normalization directly stemming from the chiral anomaly of QCD, as 
was put forward in Ref. Ill in . 

All these facts reinforce the desirability of a less model-dependent analysis of 
hadronic % decays in order to gain more insight on the hadronization of the relevant 
QCD currents. 



2. THE RESONANCE CHIRAL THEORY OF QCD 

As the energy increases approaching M p , it becomes unavoidable to explicitly introduce 
the new active degrees of freedom (the resonance mesons) into the effective action. 
The procedure, developed in Refs. |fl2L fl3 1 and known as Resonance Chiral Theory 



(RXT), is ruled by the approximate chiral symmetry of low-energy QCD for the lightest 
pseudoscalar mesons and by unitary symmetry for the resonances. Its two theoretical 
pillars are: 

• The Weinberg's Theorem [2] that ensures that writing the most general Lagrangian 
consistent with the assumed symmetries we will obtain the most general S-matrix 
elements respecting analiticity, perturbative unitarity, cluster decomposition and 
the starting symmetry principles and, in particular, local chiral symmetry in the 
presence of external fields um- 

• The large-A'c expansion of QCD. It has been suggested |fl5|1 that l/Nc (Nc is short 
for the number of colors) could serve as a perturbative expansion parameter. Many 
features of meson phenomenology find their explanation in the Nc — » 00 framework 

which supports the procedure. Relevant consequences of it are that meson 
dynamics is described, at LO in the expansion, by tree level diagrams given by 
an effective local Lagrangian including the interactions amidst an infinite number 
of zero-width resonances. However in most processes, like x decays, we need to 
include finite widths (a NLO effect in the l/Nc expansion) as we do within our 
framework II 1711 . We are also departing from the Nc — > 00 limit because we consider 
just one multiplet of resonances per set of quantum numbers (single resonance 



approximation lll8ll ) and not the infinite tower predicted. 



Once we have highlighted the key points of our study, it is time to give some details 



about the concrete way we proceed. The relevant part of the RyT Lagrangian is Ill2l.ll9l. 



20. 



+ j< + f + E + 1 *4vp E +jt/ ivppp ' (1) 

where all couplings are real, being F the pion decay constant in the chiral limit. 
The notation is that of Ref. uM- Here and in the following P stands for the lightest 
pseudoscalar mesons. Furthermore, all couplings in the second line are defined to be 



dimensionless. For the explicit form of the operators in the last line, see |20J, I2JJ, |22Q. 
In this work we have considered for the first time the VPPP vertex in the odd-intrinsic 
parity sector. Notice that we are using the antisymmetric tensor formalism to describe 
the spin one resonances. We are not considering the @(p A ) Lagrangian of Goldstone 
bosons to avoid double counting, because the LEC's of 0(p 4 ) are known to be saturated 



by resonance exchange 1130 . 

An all-important caveat, however, is that J^rxt is not QCD for arbitrary values of its 
couplings. Hence, if we want to comprehend more about the non-perturbative behaviour 
of the strong interaction, we ought to learn about the determination of the couplings 
from the underlying theory ji6\ . Guided by this purpose we will demand several known 
features of QCD to our effective Lagrangian. 

A matching procedure has been put forward [0, [L3|, [S H2L 0, [25 ] by analysing 



(three-point) Green functions that are order parameters of massless QCD and matching 
the resonance theory with the LO in their OPE expansion. Additionally, if necessary, we 
will demand to the form factors we obtain a Brodsky-Lepage-like behaviour both to the 
vector and axial- vector ones. 

Phenomenology could also help us to fix bounds on the coupling in Eq. 
For instance, Fy could be extracted from the measured r(p° — > e + e~), Gy from 
r(p° — > k + k~), Fa from T(a\ — > ny) and the A ; 's from T(a\ — > pn) which starrs on 
the T — > 3nv T processes themselves. r(co — > ny), T((o 3n) and the 0(p 6 ) correction 
to T(n — > 77) may give us information on the remaining couplings (20J]. 



3. FORM FACTORS IN x~ -> {KKn)~ v T 
T _ — > K K° n° v T can be written as 



The decay amplitudes for the processes % — > K + K it v T , x — > K°K°n v T and 



= -^Lv us u v y(l-y 5 )u^, (2) 



where 



= (p( P1 )p(p2)p(P3)\(y»-^)^ QCD \°)= (3) 

V^Ff(Q\ ShS2 )+V 2 ^(Q 2 , ShS2 ) + Q^(Q 2 , Sl , S2 ) + iV 3 ^(Q 2 , ShS2 ), 



and 



V3 l i = e t ivpaP\P2Pi » Qii = (Pi + Pi+ Pi)ii , Si=(Q-pi) 2 . (4) 

Here i*}, i = 1,2,3 correspond to the axial- vector current while drives the vector cur- 
rent. The form factors F\ and F2 have a transverse structure in the total hadron momenta, 
QV, and drive a J p = 1 + transition. The scalar form factor, F3, vanishes with the mass 
of the Goldstone bosons (chiral limit) and, accordingly, gives a tiny contribution. In dsj] 
the vector form factor did not contribute in the isospin limit, so the aim of our work is 
not only to confirm or refuse the bounds on the {A/}? 1 in that reference, but to explore 
the vector currect sector of the resonance Lagrangian. For details on the particular 



shape of the form factors, see Refs. [|21l. 12211 . 



4. ASYMPTOTIC BEHAVIOUR AND QCD CONSTRAINTS 

The explicit computation of the Feynman diagrams involved shows that the result de- 
pends only on three combinations of the four of the {c,}J =1 , two of the {di}f =1 
and three of the {g/}f =1 . The number of free parameters has been reduced from 24 to 15. 

We require the form factors of the A* 1 and currents into KKtz modes vanish at 
infinite transfer of momentum. As a result, we obtain these constraints from A^: 

F V G V =F 2 , Fy—F\ = F 2 , M 2 F 2 =M 2 F 2 , 

11 / 1F^\ 
X' = —= , ~ 0.455 , X" = - 1 5- U 7 ~ -0.094 , (5) 



2V2 f7~^ ' V f 2 



1 — — 

F 2 



contrary to X' = A, X" = and Fy = y/2F in [JsJ] . Furthermore, we have checked that 
the set of relations © are the most general ones satistying the demanded asymptotic 
behaviour in x — > 3nv T . 



Incidentally, the value of the ratio ^ can be read from [19] through: 



F 2 „ „ U a ^ an A fi-^™ lit Oh tUm.mU- F V — M A 

F 7 M^—My 



which yields -with M A = 998(49)MeV [|26|]- F v ~ 0.147 GeV. Only Xq remains free. 
It was fitted to be around 12 in Ref. ^ and derived to be some A in Ref. lfl9Tl . A fit to 
r (t — > K + K K Vx) excludes the former value and favours the second one. 

Proceeding analogously for the vector current form factor, we find that it is only 
possible to satisfy the QCD requirements provided we include the VPPP part of the 
Lagrangian in ([]])• We have found the constraints H 



The folllowing definitions have been used: C1235 = c\ + c 2 + 8C3 — C5, C125 = c\ — c 2 + C5, C1256 = 
c\ - c 2 - c 5 + 2c 6 , g 123 =gi+2g 2 -g3 and ^123 =d\ + M 2 - d 3 . 



TABLE 1. Comparison of our work with the branching ratios pub- 
lished by the experimental collaborations in units of 10~ 15 GeV. Our 
estimates for the errors are conservative educated guesses. 



Contribution/Source F(z 


-> K+K-j[-V T ) T(t~ 


-^K-K°7t Q v T ) 


Vector current, F v 


3.3 (0.4) 


2.0 (0.3) 


Axial-vector current, F A 


1.2 (0.5) 


1.0 (0.4) 


Total 


4.5 (0.6) 


3.0 (0.5) 


PDG'06 [27J 


3.465 (227) 


3.488 (454) 


BaBar'07 [28J 


3.048 (084) 




CLEO'03 [9J 


3.511 (136) 





Nr Ml , 
* = -1^2^-- ail( - - 08) ' 

o 2 = J^J^ 59 . 10 -3 (6) 

as well as C125 = 0(0) and #123 = 0, where the values in parenthesis refer to the relations 
in Ref. Il20ll . Although obtained in different kind of processes, they are quite similar. 
This points out to an error associated to the staying at LO in the large-Afc expansion 
(much) lower than the conservative 30% usually claimed. From the 24 initially unknown 
couplings in Eq. ©, only five remain free: C4, C1235, ^123, g4 and g$. 

We have also completed the computation of CO — > 3n in Ref. ||20|l by including the 



local contribution. A fit to the branching ratio 112711 fixes 2g4 ±#5 = —0.48 ±0.01, so 
only four parameters remain. In Ref. [20] the relations J123 = 0.05 and C1235 = were 
found. We have succeded in fitting the branching ratios using these values. 

These way, only C4 and a combination of #4 and #5 remain to be known. In 
T~ — > ^r _ ^ ^ v-r, the computation of the Wess-Zumino term gives zero and the VPPP 
one vanishes when imposing QCD short-distance behaviour. A fit to the branching 
ratio ITT] gives C4 = 0.06 ±0.01. This lack of two contributions translates into a 
lower branching ratio for this channel. Keeping this value, we favour g4 ~ —0.24 and 
~ gs "C g4- Remarkably, despite starting with 24 unknown couplings, we have been 
able to fix bounds on all entering the x — > KKnv x decays. 



5. RESULTS AND CONCLUSIONS 

As a result of the procedure described above, we have obtained the branching ratios 
shown in Tabled] 
Several comments are pertinent: 

• Within the very accurate SU(2) isospin limit, the following equality holds: 
r(r _ -^K + K-n-v r ) = t(t~ ^ K°K° 7t°v t y while PDG06 reports the value 

r U- ^°ZVV) = 3.624 (704) ■ 10~ 15 GeV. Taking this and the values 
presented in Table [1] into account, it is probably safer to keep a central value 



FIGURE 1. Spectral function plotted versus the momentum of the hadronic system squared for the 
channels T~ — > K K° 'n° V T (left), and T~ — > 7T~ V T (right). The dotted line corresponds to the 

axial-vector current contribution and the dashed line to the vector current one. The total contribution is 
plotted in solid line. 



close to 3.5 • 10 15 GeV but enlarge the error to some 0.5 ■ 10 15 GeV. The results 
from BaBar [|28h and CLEO Q9[| seem quite optimistic despite being rather different. 



Our results are consistent with the PDG values Il27ll . As already pointed out in 



112911 we observe a clear dominance of the vector current, which is at odds with the 



conclusions given in Refs. IIIOL BOH. 

As a new result, we find more decays to the kaon-antikaon channels. Schematically, 
we can summarize our findings as: 

r- r(tr^jr^*-v t ) 

Our prediction for the spectral functions can be seen in Fig.CQ 



This work is framed into a wider project that extends to all three meson decay channels 
of the T. We are confident that new data from CLEO, BaBar, BELLE and a future 
super-B factory might allow us to revisit and improve the works done and to help the 
completion of all others. 
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